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Aims. Four young star clusters were studied in order to characterize their anomalous extinction or variable reddening 
that could be due to a possible contamination by dense clouds or circumstellar effects. 

Methods. The extinction law (Rv) was evaluated by adopting two methods: (i) the use of theoretical expressions based 
on the colour-excess of stars with known spectral type, and (ii) the analysis of two-colour diagrams, where the slope 
of observed colours distribution is compared to the normal distribution. An algorithm to reproduce the zero age main 
sequence (ZAMS) reddened colours was developed in order to derive the average visual extinction (Ay) that provides 
the best fitting of the observational data. The structure of the clouds was evaluated by means of statistical fractal 
analysis, aiming to compare their geometric structure with the spatial distribution of the cluster members. 
Results. The cluster NGC 6530 is the only object of our sample showing anomalous extinction. In average, the other 
clusters are suffering normal extinction, but several of their members, mainly in NGC 2264, seem to have high Rv, 
probably due to circumstellar effects. The ZAMS fitting provides Ay values that are in good agreement with those 
found in the literature. The fractal analysis shows that NGC 6530 has a centrally concentrated distribution of stars 
that is different of the sub-structures found in the density distribution of the cloud projected in the Av map, suggesting 
that the original cloud has been changed with the cluster formation. On the other hand, the fractal dimension and 
the statistical parameters of Berkeley 86, NGC 2244, and NGC 2264 indicate a good cloud-cluster correlation, when 
compared to other works based on artificial distribution of points. 

Key words, open clusters: individual (NGC 2264, NGC 2244, Berkeley 86 and NGC 6530) - ISM: dust, extinction - 
Stars: pre-main sequence 



1. Introduction 

Extinction is one of the best known consequences of the 
dust present in the Interstellar Medium (ISM) . The study of 
extinction on star-forming regions reveals properties of the 
interstellar material that leads to a better understanding of 
the formation and evolution of stars. 

Each line-of-sight has an extinction law, which is char- 
acterized by the ratio of the total-to-selective extinc- 
tion Ry — Ay I E(B-V) and depends on the c ompo - 
sition of the ISM. According to ISavaee k Mathisl l| 1979ft 
its average value in the diffuse ISM is 3.1, while in 
denser regions it might reach values in the range 4 < 
Ry < 6. Anomalous high extinc tion laws have been 
found in several star-forming regions (iNeckel k Chinilll98lt 
Chini k KrugeWl983tlChim k Wargaulll990HPandev et al l 
2000: ISamal et alJl2007ft . Besides the differences of the ISM 
characteristics, high values of Ry are commonly explained 
by selective evaporation of small grains caused by the ra- 
diation from hot stars, or by grain growth in circumstellar 
environments (e.g. van den Ancker et al. 1997). The cir- 
cumstellar hypothesis seems to be supported by the cor- 
relation between reddening and extinction found by Sung 
et al. (2000) for a sample of 30 stars. Most of these ob- 
jects had their SED well fitted by assuming Ry =3.1, with 



some exceptions that show anomalous Ry and have high ex- 
tinction probably due to the circumstellar material. From 
the results of a polari metric survey i n the direction of the 
Lagoon Nebula (M8), iMcCall et all (|1990l ) estimated Ry 
= 4.64 ± 0.27 that is attributed to circumstellar effects. 
This anomalous high Ry was obtained after removing the 
normal foreground extinction. 

The dependence of the extinction law as a function of 
wavelength has also been widely discussed in the literature, 
given the interest in identifying the characteristics of dif- 
ferent types of grains causing the extinction, as well as the 
need to properly correct the reddening that affects the ob- 
servational data. In general a "universal extinction law" , 
usually a power law for wavelengths larger than 1/im, has 
been adopted. The power law (A\ oc A _/3 ) fits most ex- 
tinction curves. However, it is recognized that the index /3 
varies significantly depending on the l ine-of-sight. Based on 
data f rom the Hubble Space Telescove. iFitzpatrick k Massal 
(2009) investigated the extinction dependence on near- 
infrared (near-IR) wavelengths. They found that a "uni- 
versal" law does not apply in this case and the index /3 
has a tendency to decrease with in creasing Ry. Instead 
of ext rapolating the usual power law, IFitzpatrick k Massal 
(|2009ft suggest a different form, independent of wavelength, 
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that better describes the near-IR extinction (see discussion 
in Sect. 3.1). 

One of the most important issues directly affected by 
the correct determination of extinction is the accuracy on 
distance estimation. The cluster NGC 3603, for instance, 
has distances found in the literature varying from 6.3 to 
10.1 kpc, a discrepancy mostly due to error s on the red- 
dening correction. Adopting the normal i?y, iMelena et al.l 
420 08jj I estimated a distance of 9.1 kpc for NGC 3603, while 
a distance of 7.6 kpc is found if they assu me the anoma- 
lous e xtinction law Ry — 4.3, as proposed bv lPandev et al.l 
(2000|). The correct determination of distance is crucial in 
the estimation of physical parameters of open clusters, like 
stellar masses and ages. Particularly in star-forming regions 
special attention must be paid to the probable occurrence 
of anomalous or differential extinction. 

Several methods can be employed to estimate the in- 
terstellar extinction in the direction of young star clusters, 
which can be due to different effects: (i) the foreground in a 
given line-of-sight that is pervaded by interstellar material 
in between the cluster and the observer; (ii) the presence of 
a dark cloud associated to the cluster, (iii) the individual 
extinction caused by circumstellar material. 

The use of U-B x B-V colour-colour diagram is the clas- 
sical method for estimating the average extinction towards 
open clusters (e.g. Prisinzano et al. 2011) and it is specially 
useful on the lack of spectroscopic observations (e.g. Jose 
et al. 2011). 

In addition to colour-colour diagrams in the UBV 
bands, two-colour diagrams (TCDs) are also interesting for 
the study of extinction as proposed by Chini & Wargaw 
(1990). These diagrams are presented in the form V-A x 
B-V, where A refers to different bands. The distribution of 
the stars in TCDs is roughly linear. Anomalies in the ex- 
tinction law are determined from the comparison between 
the distribution of field stars (which follows a normal ex- 
tinction law) with the distribution of the cluster members 
(w hich may follow a n anomalous extinction law). 

iJose et ail (|2008f ) adopted TCDs to verify the extinction 
law in the direction of the cluster Stock 8, for instance. They 
found the same slope for the inner region of the cluster (r < 
6') as well as for a larger area (r < 12') that also includes 
field s tars. Another examp le is the cluster NGC 3603 for 
which IPandev et a l. (2000) estimated Ry = 4.3, indicating 
a remarkable difference between the distribution of field 
stars and members of the cluster on the TCDs. 

The aim of the present work is to characterize the ex- 
tinction in the direction of different young star clusters, 
by determining the extinction law and searching for possi- 
ble spatial variations. We selected four well-known clusters, 
whose characteristics are compiled in the Handbook of Star 
forming regions edited by Bo Rcipurth (2008): B erkeley 86 
is fou nd in the Cygnus OBI region (|Reipurth fc Schneider! 
20081): NGC 2244 is associa ted to the Rosette Nebula 

NGC 2264 is related to the 
and NGC 6530 is lo- 



Table 1. List of the studied clusters. 
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cated near to the Lagoon Nebula |Tothill et al.l l2008). 

The motivation in choosing these well-known objects is 
to refine the previous extinction determinations, by adopt- 
ing same criteria for selection and analysis of data sets, 
in order to compare our results with the characteristics of 
the clusters environments. The paper is organized as fol- 
lows. Section 2 presents the sample and the information 
available in the literature for the members of the clusters 
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and the clouds associated to them. Different methods are 
adopted in Sect. 3, aiming to estimate Ry in the direction 
of the clusters. Section 4 describes an automatic method 
to fit the ZAMS reddened colours to the observed data, 
providing an accurate estimation of visual extinction. In 
Sect. 5 we develop an analysis of the fractal dimension of 
the clouds by comparing the spatial distribution of cluster 
members with statistical parameters related to clustering. 
The discussion of the results and the conclusions are pre- 
sented in Sect. 6. The colour diagrams utilized for studying 
the extinction are presented in the Appendix A. 



2. Description of the Sample 

The list of clusters and their main characteristics are given 
in Table 1 . This section is dedicated to summarize the infor- 
mation found in the literature, and to describe the adopted 
criteria in selecting the cluster members and their available 
observational data. We also performed an analysis of visual 
extinction maps and far-IR images of clouds against which 
the clusters are projected. 



2.1. Selected Clusters 

Berkeley 86 is a particularly small cluster associated to 
the Cygnus OBI region. Its youth w as revealed by th e 
presence of O typ e stars, discovered b y [Sanduleakl (|1974[ ). 
Recent results by iBhavva et al.l (|2007l ) indicate a distance 
of 1585±160 pc. They suggest that a low level star forma- 
tion episode occurred 5 Myr a go, and another more vigorous 
started during the last 1 Mvr. lReipurth fc Schneider] (2008) 
presents a comparison of gas and dust distributions in the 
Cygnus X region, which is influenced by the UV radiation 
from OBI association (see their Fig. 3). They show that 
Berkeley 86 is isolated, projected against an area totally 
free of gas and dust, suggesting low levels of extinction for 
this cluster. Yadav and Sagar (2001) argue that Berkeley 86 
suffers non-uniform reddening with colour-excess varying in 
the range 0.24 < E(B-V) < 1.01, while Forbes (1981) sug- 
gests a more uniform estimation of E(B-V)=0.96±0.07. 

NGC 2244 is located at the centre of the Rosette Nebula. 
It is a prominent OB association probably responsible for 
the evacuation of the central part of the nebula. Th e first 
photometric study of this cluster (jJohnsonl Il962f) indi- 
cated a colour-excess E (B-V) = 0.46 a ssum ing Ry=3.0 that 
was la ter confirmed by iTurneil (|1976ft and lOgura fc Ishidal 
(198l|). Perez et al. (1987) found anomalous Ry for some of 
the cluster members and suggested the coexistence of main 
sequence (MS) and pre-MS stars. Roman- Zun i ga and Lada 
(2008) suggest ages of 3±1 Myr. IWang et al.l (|2008h stud- 
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ied the X-ray sources detected by Chandra in the Rosette 
region. They verified a double structure in the stellar radial 
density profile, suggesting that NGC 2244 is not in dynam- 
ical equilibrium. 

NGC 2264 is one of the richest cluste rs in terms of mass 
range, presenting a well defined pre-MS (Dahm 2008). The 
estimated ages vary from 0.1 to 10 Myr (Flaccomio et al. 
1997, Rebull et al. 2002). Surveys in Ha and X-rays re- 
vealed the presence of about 1000 members. The cluster 
is observed projected against a molecular cloud complex 
of ~ 2° x 2° and is located at 40' north of the Cone 
Nebula. The interstellar reddening on thi s cluster is be- 
lieved to be low: ,4y=0.25 estimated by IWalkerl (|1954j ) 
ad opting E(B-V)= 0.08 2 and i?y=3.08. A more recent work 
bv iRebull et all (|2002f ) derived ^4y=0.45 fr om a spectral 
study in a sample with more than 400 stars. lAlencar et al.l 
(2010) used the CoRoT satellite to perform synoptical anal- 
ysis of NGC 2264 with high photometric accuracy. Their re- 
sults suggest dynamical star-disc interaction for the cluster 
members, indicating they are youn g accre ting stars. 

In the study by iTothill et all (|2008h of the Lagoon 
Nebula and its surroundings, the main emphasis is given to 
NGC 6530, a star cluster associated to an HII region that 
lies at about 1300±100 pc. The estimated age of NGC 6530 
is less than 3 Myr (Arias et al. 2007). Despite being in 
a line-of-sight that contains a high concentration of gas, 
NGC 6530 appears decoupled from the molecular cloud 
once the cluster members do not seem to be very embed- 
ded. This is consistent with evidenc es that NGC 6530 i s 
inside the cavity of the HII region (jMcCall et al.l Il990h . 
indicated by measurements of the expanding gas (Welsh 
1983). A range of 0.17 to 0.33 mag has been reported for 
the colour-exc ess in the direction of NGC 6530 (Mayne & 
Naylor, 2008). Ivan den Ancker et alJ (|l997h found a normal 
extinction law for the intra-cluster region, while Arias et al. 
(2007) suggest Ry = 4.6 for some of the embedded stars. 
Anomalous extinction has been also reported for cluster 
members individually, for instance the star HD 164740, for 
which Fitzpatrick & Massa (2009) suggest E(B-V) = 0.86 
and 5.2 < R v < 6.1. 

2.2. Data extraction and Cluster membership 

Two main catalogues of open clusters can be found in the 
literature: (i) the WEBDAEI database and (ii) the cata- 
logue DAML023 by iDias et all (|2002bh with the compila- 
tions: Tables of members h ip and mean proper motions for 
open clust e rs, IDias et alJ (|2002aft based on TYCH02 and 
IDias et alJ (|2006l) based on the UCAC2 catalogue. 

From DAML02 we searched for clusters with ages up 
to 10 Myr and distances up to 2 kpc, only selecting those 
which had members with UBVI photometry and spectral 
classification available in the WEBDA database. A sum- 
mary of the characteristics of the sample selected from these 
two databases is listed in Table [1] 

In order to search for relevant information of stars 
located in the direction of the clusters, three data sets 
were extracted from WEBDA: (i) UB photometry, (ii) VI 
photometry, and (iii) equatorial coordinates (J2000). Each 
source has the same identification in all of these data sets, 
but the photometry is originated from different works. For 

1 www.univie.ac.at/webda/navigation.html 

2 www.astro.iag.usp.br/~wilton 
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Fig. 1. Normalized number of stars distributed according 
to membership probability (P%), which is based on their 
proper motion. A dotted line indicates the minimum prob- 
ability (Pmin) that was adopted to separate the most prob- 
able members from field stars. Top: Berkeley 86 (left) and 
NGC 2244 (right). Bottom: NGC 2264 (left) e NGC 6530 
(right). 



the consistency in the data analysis, we selected only clus- 
ter memberships having UBVI photometry provided by a 
single work. The number of studied stars and the corre- 
sponding references to the photometric catalogues are listed 
in Table 2, which also gives information on the number of 
stars having avai lable spectral data. 

According to lYadav fc Sagar (2001) the correct identi- 
fication of cluster members is crucial for the assessment of 
extinction in that direction and the most reliable selection 
is based on kinematic studies (proper motion and radial 
velocity). DAML02 provides a list of stars with JHK pho- 
tometry - 2MASS catalogue (jCutri et alJl2003f ) - and also 
the probability of the star belonging to the cluster (P%), 
which is determined from its proper motion. In order to 
verify if P% is available for our WEBDA sample, their co- 
ordinates were compared with those listed by DAML02. A 
coincidence of sources was only accepted for objects having 
less than 5 arcsec of difference between coordinates. 

The distribution of the number of stars as a function of 
P% was evaluated on basis of histograms that are presented 
in Fig. 1. A bimodal distribution is verified, enabling us to 
distinguish between member s and possible field stars. As 
proposed bv lYadav fc Saeail (|200lh . the pollution by field 
stars should be reduced if only stars with high values of P% 
are considered members of the cluster. Based on Fig. 1, for 
each cluster a minimum probability (P m m) was adopted 
to separate field stars, as indicated by dotted lines in the 
histograms. 

The spatial distribution of the stars was also checked 
looking for possible preferential concentrations, as a func- 
tion of P%. However, no trend was found in the members 
distribution. 

2.3. Comparison with Molecular Clouds 

Extinction caused by interstellar dust found in the line- 
of-sight can be determined by models that reproduce the 
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Fig. 2. Map of visual extinction over-imposed by contours measured from IRIS-IRAS images at 100/im band. The 
contours vary from 50 - 350 MJy/Sr with steps of 100 MJy/Sr, excepting for NGC 6530 that has contours starting at 
500MJy/Sr with steps of 1000 MJy/Sr. The area of the cluster is indicated by the black central square. 



stars counts in the Galaxy (e.g. Amores & Lepine, 2005) 
or maps of visual extinction (Gregorio-Hetem et al. 1988; 
Schlegel et al. 1998; Dobashi et al. 2005, among others). 
More recently, colour-excess maps obtained from near-IR 
surveys have been used to derive the visual extinction in 
the direction of molecular clouds, as for instance the Trifid 
nebula, studied by Cambresy et al. (2011) based on 2MASS, 
UKIDSS, and Spitzer data. 

Figure 2 displays the position of each cluster projected 
against the visual extinction map obtained from the 2MASS 
data (K band) . The Catalogue of Dark Clouds by Dobashi 
et al. (2005) was used to extract the maps for the regions 
containing our clusters. In order to compare the extinction 
with the far-IR flux density of the clouds, the Ay maps are 
over-imposed by contours of IRAS-IRIS data at lOO^rn. A 
significant variation of both, Ay and far-IR emission, can be 
found in the direction of the clusters, excepting Berkeley 86 
that appears projected against a uniform field. High levels 
of far-IR emission are found in the direction of the clusters 
NGC 6530 and NGC 2264, which coincide with dense re- 
gions in the Ay maps. Table 1 gives the ranges of flux at 
100/im and Ay in the direction of the clusters. In spite of 
the sub-structures of the clouds being more evident in K 
band, we adopted the optical (DSS) maps from Dobashi et 
al. (2005) to measure Ay. 

3. Study of the Extinction Law 

3.1. Ry estimation based on colour-excess 

One of the methods to obtain Ry uses colour-excess, 
which can be determined for each cluster member that 
has well-known spectral type. Intrinsic colours were 
adopted from Bessel et al. (1998), based on surface 
gravity and effective temperature that were r espec- 
tively extracted from Straizvs &; Kurilie"nel (|1981l ) and 
Ide Jager fc NieuwenhuiizenT ! 19871) as a function of the 
spectral type and luminosity class available for some stars 
of the sample. 

The expression -j^ = a(x) + was adopted from 



iCardelli et al.l (|1989f ) in order to obtain relations between 
Ry and colour-excess. The results are similar to those based 
on van de Hu lst's theoretical extinc t ion (e .g. Johnson 1968). 
According to iFitzpatrick fc Massal (|2009l ). these theoretical 
curves seem to be valid only for low- to moderate values of 



Ry , which tend to be underestimated in the case of anoma- 
lous high extinction laws. They suggested a new formula 
Ry[K] = 1.36 ^|g~yj — 0.79, which better reproduces the 
extinction law in different lines of sight. 

Due to the lack of spectral type information for our 
entire sample, we did not determine Ry for each cluster 
member. Instead of that, we examine the expected ranges 
of Ry as a function of the distribution of colour-excesses 
compared to the theoretical extinction laws on the E(V-K) 
x E(B-V) plot presented in Fig. 3. An illustration of the 
dispersion found for Ry [K]=6 is shown by th e hatched area, 
corres ponding to the extinction laws given bv lCardelli et all 
(|1989l ) (upper line) and by Fitzpatrick & Massa (2009) 
(lower line). For Ry [K]<4 the theoretical lines from both 
works are quite similar. 

Considering the few objects having known spectral type, 
this analysis cannot be conclusive for our sample, but gives 
some indication on the variation of E(B-V) and Ry for each 
cluster. NGC 2264 shows the largest variation, with E(B-V) 
ranging from 0.01 to 0.3 mag as measured for a sub-sample 
of 20 stars. A degenerate distribution of Ry is found in 
this case, spreading from 2.1 to > 6. It can be seen in Fig. 
3 that NGC 2264 seems to present a bimodal distribution, 
where part of the members follows the normal extinction 
law, while another part has Ry > 6. 

There are several possible causes to the large variations 
found in NGC 2264: (i) observational problems that would 
give inaccurate optical photometry or spectral type deter- 
mination; (ii) actual variations in each line-of-sight around 
the cluster area; (iii) anomalies caused by circumstcllar 
environment. In principle, the two first options may be 
disregarded, since the observational information has been 
checked in other data sets. Furthermore, there are no corre- 
lation between the position of a given cluster member and 
its Ry . We conclude that a presence of circumstellar matter 
seems to better explain the anomalies. It can also be noted 
in Fig. 3 that E(V-K) is particularly large for seven objects 
of NGC 2264, inconsistent with the low E(B-V) values, sug- 
gesting K band excess that is expected in the presence of 
circumstellar matter. Six of these stars have spectral type 
later than AO, for which it is important to take into account 
the errors of spectroscopic data that may cause significant 
uncertainties on colour excess (see Sect. 3.3). However, it is 
interesting to note that some of these objects also show an 
offset from the normal distribution of the observed colours 
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on TCDs (see Sect. 3.2 and Fig. A2), which is more consis- 
tent with individual anomalies of Ry than possible errors 
on photometry. 

Smaller variations are found for the other clusters: mem- 
bers of NGC 2244 have R V [K] ~ 3.1; NGC 6530 shows 
a trend to slightly larger values (3.1< i?y[K]<4), while 
Berkeley 86 tends to have lower values (2.7<Ry[K]<3). 
However, these results are based on few stars and need to 
be confirmed by different methods of Ry estimation. 



3.2. Two-Colour Diagrams 

In order to search for anomalies in the extinction law, 
we also constructed TCDs following the method proposed 
by Chini et a l . (198 3 . 1990) and more recently used by 
iPandev et all (|2000[) . jjose et al.l (|2008D . Chauhan et al. 
(2011) and Eswaraiah et al. (2011), among others. TCDs 
can be used to perform a qualitative analysis of the nature 
of the extinction by using the relation: R c = R n x (a c /o„), 
where and R n is the usual normal extinction law, a c is the 
slope of the linear fit for the cluster members, and a n is the 
slope for field stars. 

TCDs use V-A colours, where A represents the IJHK 
bands. In these diagrams the intrinsic colours show a linear 
distribution up to B-V ~ 1. This limit defines a restricted 
spectral range of validity, where the ZAMS distribution 
corresponds to a straight line. In the TCDs, the distribu- 
tion of field stars is well reproduced by the ZAMS intrinsic 
colours, once they are reddened by using the mean E(B-V) 
estimated in the direction of the cluster. In this case, the 
ZAMS fitting is obtained with i?y=3.1, which we adopted 
to simulate the distribution of the field stars that must be 
compared to the cluster members. 

The fitting of our sample is restricted to the cluster 
members having P% > P m im aiming to avoid confusion 
with field stars. Since the TCD analysis is only valid for 
a very narrow range of B-V, our calculations are based on 
objects with spectral type earlier than B8. For each studied 
band, the slope of the linear fit for the selected cluster mem- 
bers (a c ) was compared to the slope of the normal ZAMS 
distribution (a„). 

Figures Al and A2 show the diagrams V-I, V-J, V-H, V- 
K x B-V, for which a„ is 1.09±0.02, 1.87±0.03, 2.42±0.03 
and 2.51 ±0.04, respectively. The effect of reddening com- 
bined with anomalous extinction is illustrated in the V-K 
x B-V diagrams, where extinction vectors were plotted to 
represent Ay = 2 mag for two choices of Ry, as used by 
Da Rio et al. (2010), for instance. 

Table [2] presents the results from the TCD analysis giv- 
ing the fitted line slope in each band, which is used to 
estimate the mean Ry obtained for the early type stars. 
Within the estimated errors, a good agreement is found 
in all bands, excepting NGC 6530. In this case Ry ranges 
from 4.5 to 6.2. This is the only cluster of the sample clearly 
showing anomalous extinction law with significant depen- 
dence on wavelength. Berkeley 86 and NGC 2264 seems 
to show a normal law. As verified in Sect. 3.1 (Fig. 3), 
NGC 2244 also has a bimodal distribution in the TCDs: 
part of the objects follows the normal law, while other part 
shows anomalous extinction. 




E(B-V) 

Fig. 3. Different extinction laws compared to the colour- 
excess distribution for cluster members that have available 
spectral type. Error bars illustrate the maximum uncer- 
tainty expected, if the spectral type is not well determined, 
for a G8 star. 



3.3. Errors estimation 

Since our work is based on photometric and spectroscopic 
data from literature, we have evaluated the global uncer- 
tainties taking into account the different sources of errors. 

According to the references that provided the photo- 
metric data (see Tabled]), the error on the magnitudes and 
colours is ^0.08 mag in the worst case. Since this value is 
smaller than the symbol size on the colour-colour diagrams 
(see Figs. Al and A2), it does not not affect the linear fit- 
ting used to estimate the variations on the extinction law. 

Most objects with known spectral types in our sample 
are O or B type stars, whose errors in the characterization 
of spectral types are not significant in our calculation of 
colour-excess. The only exception is NGC 2264 that has 
some A, F or G stars for which we adopted a maximum 
uncertainty of two spectral subtypes. 

By combining both, spectroscopic and photometric er- 
rors, in order to derive colour-excess, we estimate the max- 
imum deviation (the worst case) of 0.3 mag on E(B-V) 
and 0.4 mag on E(V-K) for a G8 type star of NGC 2264. 
Considering that we do not use the colour-excess to derive 
the extinction law, only performing a qualitative evaluation 
of the expected range of Ry, these errors do not affect the 
discussion based on colour-excess. 



4. Main Sequence Fitting Algorithm 

The fit of photometric data to the theoretical isochrones 
is commonly used to determine star clusters parameters. 
In order to avoid subjective results that could be derived 
from visual inspection, different methods of automatic fit- 
ting have been p roposed to optimize th e solution for multi- 
variable fittings. [Monteiro et al.l (|2010f l. for instance, devel- 
oped a technique based on the cross-entropy global opti- 
mization algorithm. This method proved to be a powerful 
tool to fit the observed data in c o lour-m agnitude diagrams. 
More recently, iMonteiro fc Diasl (|201lf l used this technique 
to determine open-cluster parameters by using BVRI pho- 
tometry, which provides reliable estimation of distance, ex- 
tinction, mass and age. 
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Table 2. Extinction laws from literature and results obtained from the TCD analysis. 



Cluster 


N tot 




Kv 


OellJ 


Rvll] 


a c [,T] 


Rv [JJ 


a c [Hj 


RWH] 


Qc[K] 


R V [KJ 


Berkeley 86 


16" 


3" 


3.0* 


1.20 


3.39 ± 0.11 


1.98 


3.29 ± 0.11 


2.24 


2.87 ± 0.10 


2.41 


2.98 ± 0.10 


NGC 2244 


288 b 


19 f 


3.2 - 3A j 


1.14 


3.22 ± 0.11 


2.02 


3.36 ± 0.13 


2.49 


3.19 ± 0.10 


2.54 


3.13 ± 0.22 


NGC 2264 


179 c 


24 s 


3.1 - 5.2 fc 


1.04 


2.94 ± 0.20 


1.91 


3.17 ± 0.08 


2.39 


3.07 ± 0.08 


2.50 


3.09 ± 0.11 


NGC 6530 


24 d 


8 h 


3.1 - 5.4 ! 


1.59 


4.49 ± 0.14 


3.45 


5.73 ± 0.38 


4.31 


5.52 ± 0.19 


5.01 


6.19 ± 0.32 



Notes: Ntot is the total number of studied stars and Nst refers to those having available spectral type, the respective references to photometric and 
spectroscopic catalogues are: (a) Deeg & Ninkov (1996); (b) Park & Sung (2002); (c) Sung et al. (1997); (d) Sung et al. (2000); (e) Forbes et al. (1992); 
(f) Johnson (1962), Morgan et al. (1965), Conti & Leep (1974), Hoag & Smith (1959), Wolf et al. (2007); (g) Morgan et al. (1965), Walker (1956), 
Young (1978); (h) Hiltner (1965). The extinction law references are: (i) Bhavya et al. (2007); (j) Ogura & Ishida (1981), Perez et al. (1987); (k) Walker 
(1956), Perez et al. (1987); (1) Neckel & Chini (1981), Arias et al. (2006). 



Since we are interested in solving a single variable prob- 
lem, we decided to adopt a simple method of fitting. In 
order to improve the efficiency in determining the visual 
extinction for the clusters, we developed an algorithm that 
automatically fits the theoretical ZAMS curve to the obser- 
vational data in a colour-colour diagram. The best fitting 
provides the average Ay in the direction of the clusters. 

We define Zq \ubj , bvj] as the curve formed by the unred- 
dened ZAMS (jSiess et al.ll200rl on the U-B x B-V plan (CC 
diagram). The ZAMS is defined by i=l..n points (intrinsic 
colours), while the distribution of observed colours is ex- 
pressed by (J, where crj[xj,yj} is the position occupied by 
the star j on the CC diagram, with j—l..m. In this case, m 
is the total number of cluster members. 

In order to obtain a reddened ZAMS curve Za v [ub, bv] 
by adding the extinction Ay to Zo[u6i,6fi]. The extinc- 
tion was calculated by using the relation ^ adopted from 
ICardelli et all t)l989f ) and the Ry determined from the TCD 
analysis (see Sect. 3). The reddened ZAMS was obtained 
by using R v = 3.1 for NGC 2264, NGC 2244 and Berkeley 
86. In the case of NGC 6530 we adopted Ry = 4.5 that was 
estimated from the V-I x B-V diagram. 

Following the method proposed by Press et al. (1995), a 
Cubic Spline interpolation is applied to the reddened ZAMS 
in order to provide a theoretical curve to be fitted to the 
observational data. The best fitting is achieved by searching 
for an Ay that minimizes the distance modulus between the 
data and the interpolated ZAMS curve, given by: 



E 



{ubj - Xj f + (bvj - y.jf 



where [ubj,bvj] represents the Za v point that is closer to 

Figure A3 presents the fitting of the observed colours 
compared to the curve of intrinsic colours on the U-B x 
B-V diagram. It must be kept in mind that Ay obtained 
from the ZAMS fitting is only an average. A clear dispersion 
around this average can be noted in Fi g. A3 for NGC 2264 
and NGC 2244. According to iBuTkl (fl975h . a maximum 
dispersion of AE(B-V) = 0.11 may be due to other effects 
like duplicity, rotation and age differences. This expected 
dispersion is illustrated by two curves displayed in the bot- 
tom panels of Fig. A3. In this case, only cluster members 
with P%>P m m are plotted. 

Clusters with uniform or non- variable reddening present 
a distribution of observed colours falling into the expected 
range of dispersion. The same cannot be said for several 
members of NGC 2264 and NGC 2244 that have colours out 
of AE(B-V) = 0.11, which suggest non-uniform reddening. 



In summary, the results from the study of extinction law 
and the main sequence fitting indicate that Berkeley 86 
seems to suffer uniform reddening and normal extinction 
law. NGC 2264 and NGC 2244 also present normal extinc- 
tion in average, however the reddening is variable for several 
members of both clusters. NGC 6530 shows anomalous ex- 
tinction (high Ry) and uniform reddening, as indicated by 
the small variation of E(B-V). 

Our results on reddening are in good agreement with 
those from literature, as can be seen in Table 3 that gives 
the colour excesses estimated by: (i) determining the mean 
value of E(B-V) on basis of spectral type, as described in 
Sect. 3.1, and (ii) the ZAMS fitting. 

In order to verify if the colour-excess measured for the 
stars is compatible with the extinction effects due to the 
cloud present in the same line-of-sight, we also present in 
Table 3 an estimation of E(B-V) that was converted from 
the extinction map by adopting Ry = 3.1 and assuming 
the minimum and the maximum Ay values given in Table 
1. 

Since all objects of our sample are very young, they 
still are physically associated to the cloud complexes where 
they have been formed, but none of them appears to be 
deeply embedded. The low extinction of the clouds (Ay < 
2.5 mag) in the direction of Berkeley 86 and NGC 2264 (see 
Table 1) gives E(B-V) similar to the results found in the 
literature (respectively Forbes 1981 and Johnson 1962), as 
well as those obtained by us, confirming that they are not 
surrounded by dense material. When comparing NGC 2264 
and NGC 6530 with their respective clouds, it can be noted 
that both are projected against regions showing high levels 
of dust emission. The flux at 100 /im reaches 720 MJy/Sr 
in the case of NGC 2264 and 5000 MJy/Sr for NGC 6530, 
corresponding to high extinction (Ay = 3.8 and 4.2 mag, 
respectively) that suggest the presence of large amounts 
of cloud material. However, both clusters are foreground 
to these regions, as indicated by the low level of interstel- 
lar reddening previously estimated for NGC 2264 (Walker 
1956, Perez et al. 1997, Sung et al. 1997, Rebull et al. 2002) 
and NGC 6530 (Mayne & Taylor 2008) and confirmed in 
the present work. 

Complementing the study of interstellar reddening, we 
also performed a fractal analysis that compares the param- 
eters of the clouds with the spatial distribution of cluster 
members, which has been suggested to be a quantitative 
method to discuss the relation of environmental conditions 
with the origin of star clusters, as described in the next 
section. 
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Table 3. Colour-excess E(B-V) estimated for the clusters. 





Berkeley 86 


NGC 2244 


NGC 2264 


NGC 6530 


Literature 


0.96±0.07 a 


0.46" 


0.06 - 0.15 c 


0.17 - 0.33 d 


S.T. 


0.98±0.04 


0.42±0.06 


0.13±0.03 


0.33±0.02 


ZAMS 


0.89±0.06 


0.44±0.06 


0.05±0.06 


0.21±0.06 


A v Map 


0.55 - 0.81 


0.12 - 0.71 


0.28 - 1.22 


1.03 - 1.35 



Description of Column 1: Literature: E(B-V) given by (a) Forbes (1981), 
(b) Johnson (1962), (c) Perez ct al. (1987) and Rebull et al. (2002), (d) 
Maync & Taylor (2008); S.T.: mean results based on spectral type (see 
Sect. 3.1); ZAMS: main sequence fitting (see Sect. 4); Av Map: visual 
extinction, from Dobashi ct al. (2005) map, converted into E(B-V) by 
adopting normal Rv ■ 

5. Fractal statistic 

Aiming to investigate a possible correlation between the 
projected spatial density of the cluster and its correspond- 
ing cloud, we perfor med a fractal analys i s base d on the tech- 
niques presented by iHetem fc Lepind jj-993). Our results 
we re compared to art ificial clouds and clusters simulated 
bv lLomax et al.l ()201ll ) (hereafter LWC11). Both works dis- 
cuss the behaviour and evolution of molecular clouds under 
the interpretation of fractal statistic. The fractal dimension 
measured on realistic simulations of density structures is 
expected to be related to physical parameters, which are 
density dependent like cooling function, dissipation of tur- 
bulent energy and Jeans limit. However, despite the statis- 
tical similarity between fractal and actual clouds, the link 
between geometry and physics still relies on empirical con- 
cepts. 

5.1. The perimeter-area relation of the clouds 

The clouds studied by LWC11 are artificial, with density 
profile given by: 



p(r) = Po 




where po is the density at r = tq, and a defines the den- 
sity law. A possible interpretation of a is related to th e 
definition of fractal dimension given bv lMandelbrot](|l983ll : 
NM{r) = [^] Dm , where N m(t) represents the number of 
self-similar structures observed at scale r < tq. 

LWC11 adopted a definition of fractal dimension similar 
to the capacity dimension given by D = l -j§^jj^r, where 

N(r) is the number of regions of effective side r occupied 
by data points. A set of points has fractal characteristic 
if the relation logN(r) x log(l/r) tends to be linear in a 
given range of r. The capacity dime nsion is determ ined by 
the slope of this linear distribution (|Turcottdll997l ). 

Our statistical analysis of the clouds is based on the vi- 
sual extinction maps shown in Fig. 2. The fractal dimension 
(D2) of contour levels is me asured by using t h e peri meter- 
area method described by IHetem fc Lepind (|1993T ): p oc 
aP" 2 ! 2 , where p is the perimeter of a given Ay contour 
level and a is the area inside it. 

The perimeter-area dimension depends on the resolu- 
tion of the maps and their signal to noise relation, (Sanchez 
ct al. 2005). We studied 2°x 2° regions covered by 123x123 
pixels, which gives to all the maps the same resolution (~ 
0.02° per pixel). Since NGC 2264 is less distant than the 
other clusters, by a factor of about 2, we have degraded 



the adopted Ay map in order to simulate a distance of d~ 
1500 pc, which gives a small difference in the measured frac- 
tal dimension (within the error bar). In order to minimize 
errors due to low signal-to- noise ratio (S/N), the lower con- 
tour level adopted in the calculation of D 2 was chosen to 
provide S/N > 10. By this way, the low density regions of 
the maps were avoided in the calculations. 

5.2. The Q parameter of the clusters 

Since stars are formed from dense cores inside molecular 
clouds, the spatial distribution of young cluster members is 
expected to be correlated to the distribution of clumps. This 
correlation can be inferred from the fractal dimension mea- 
sured on the cloud compared to the Q parameter measured 
on the cluster. In the technique proposed by Cartwright & 
Whitworth (2004), Q is related to the geometric structure of 
points distribution and statistically quantifies fractal sub- 
structures. Studies on the hierarchical structure in young 
clusters have used the Q parameter to distinguish frag- 
mented from smooth distributions (e.g. Elmegreen 2010). 

Two parameters are involved in the Q estimation: m, 
the mean edge length, which is related to the surface den- 
sity of the points distribution, and s that is the mean sep- 
aration of the points. Distributions with large-scale radial 
clustering, which causes more variation on s than to, are 
expected to have Q > 0.8. On the other hand, Q < 0.8 
indicates small-scale fractal sub-clustering, where the vari- 
ation of to is more important than s. 

The dimensionless measure Q is given by Q = m/s, 
with 

and 

2 N 
* = N(N-1)R N H l^^l 

v ' JV i=l j=l+l 

where N is the number of points in the set, rrii is the edge 
length of the minimum spanning tree, and is the po- 
sition of point i. The area A jv corresponds to the small- 
est circle encompassing all points, with radius defined by 
Rn = (4r L ) 5 ' To construct t he minimum sp anning tree, we 
used the algorithm given by iKruskall (|1956i ) and to deter- 
mine the smallest circle e ncompassi n g all p oints we adopted 
the method proposed bv iMegiddol (|l983h . This technique 
was adopted by us and employed to the spatial distribu- 
tion of the cluster members. 

The uncertainties in the estimation of to, s and Q were 
calculated by using the bootstrapping method presented by 
Press et al. (1995). This technique uses the actual set of po- 
sitions of the cluster members - Sq with N data points - to 
generate a number M of synthetic data sets - Si, 6>2 . . -Sm 
- also having N data points. A fraction / = 1/e ~ 37% of 
the original points is replaced by random points within the 
limits of the original cluster. For each new set, the parame- 
ters to, s and Q were calculated by adopting M=200. From 
these measurements, we derived 1 a deviations Ato, As and 
AQ that are presented in Table 4. 
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Table 4. Statistical parameters obtained from fractal anal- 
ysis 



Cluster 


Q 


m 


s 


D 2 


Berkeley 86 
NGC 2244 
NGC 2264 
NGC 6530 


0.73±0.12 
0.75±0.02 
0.76±0.02 
0.85±0.11 


0.69±0.12 
0.59±0.02 
0.61±0.05 
0.60±0.10 


0.93±0.14 
0.78±0.02 
0.81±0.07 
0.70±0.11 


1.39 ± 0.04 
1.37 ± 0.05 
1.47 ± 0.04 
1.34 ± 0.03 







□ 


2264 


CO 




A 


6530 






X 


//-.-. '. 






O 





5.3. Comparing clouds and clusters 

Considering that our calculations are made over two dimen- 
sion maps, while LWC11 used projections of 3D images, the 
comparison of measurements of fractal dimension can be 
done by adopting the equivalence D2 ~ D - 1. 

Figure 4 shows the Q parameter as a function of fractal 
dimension obtained by LWC11. In order to illustrate the 
offset of Q due to differences on resolution (or number of 
points) , the clustering statistics for fractal distributions (D 
— 2.0 - 3.0) of artificial data is displayed for two data sets: 
N=1024 points and N=65536 points. For comparison with 
our results, the hatched area in Fig. 4 represents the error 
bars corresponding to the lower resolution data set. 

Despite the low number of members studied in our clus- 
ters, their distribution are consistent with the LWC11 set 
of N=1024 points, excepting NGC 6530 that tends to be 
displaced from the cloud-cluster expected correlation. 

Cartwright & Whitworth (2004) suggested that clusters 
showing central concentrated distributions have Q increas- 
ing from 0.8 to 1.5. On the other hand, clusters with fractal 
sub-structures have Q decreasing from 0.8 to 0.45 related 
to fractal dimension decreasing from D2 = 2 to D2 = 0.5 
(we are using here the same relation between 2D images 
and 3D distributions adopted above). In fact, NGC 6530 
is the only cluster of our sample showing indication of ra- 
dial profile distribution of stars (Q = 0.86), which is not 
compatible with the fractal structure of its projected cloud 
(Z?2 ~ 1-3). The other clusters show members distribution 
with fractal sub-clustering structure, as indicated by their 
Q < 0.8 that is in agreement with the fractal dimension 
measured in the respective projected clouds. 

Several works have statistically proved the correlation 
between fractal dimension, estimated from cloud maps, 
and Q parameter, measured on cluster stars distributions. 
Camargo et al. (2011), for instance, summarizes the con- 
cepts first discussed by Lada & Lada (2003) suggesting that 
the structure of embedded clusters is related to the struc- 
ture of their original molecular cloud. Fractal structures are 
observed in clouds showing multiple peaks in their density 
profile (Cartwright & Whitworth 2004, Schmeja et al. 2008, 
Sanchez et al. 2010, LWC11). 

Based on the comparison with artificial data, we suggest 
for our clusters a correlation between the fractal statistics of 
the cloud and the cluster members distribution, excepting 
for NGC 6530. In this case, the measured fractal dimension 
does not indicate an uniform density distribution in the 
cloud, which is expected for associated clusters having cen- 
tral concentrated distributions. For this reason, our fractal 
analysis indicates that NGC 6530 had a different scenario 
of formation, when compared to the other three clusters. 
Our argument for discussing the cluster formation comes 
from the interpretation of the cloud-cluster relation, which 
is derived from the fractal analysis. The meaning of this 
relation is the comparison of the physical structure of the 




Fig. 4. Comparison of our clusters with the Q parameter 
and fractal dimension obtained by LWC11 for artificial data 
by using different sets of points: N=65536 (dashed line) and 
N=1024 (full line and respective deviations shown by the 
hatched area). Symbols with error bars are used to show our 
results that were estimated from the spatial distribution of 
cluster members (Q) and projected clouds (-D2). 



cluster with its remnant progenitor cloud. By this way, the 
cloud-cluster relation may tell us about the original gas dis- 
tribution of the cloud that formed the cluster. It also may 
give us information on how that particular cloud probably 
evolved since forming the cluster. 

The concentrated distribution of NGC 6530 indicates an 
original material more concentrated than the fractal struc- 
ture of the cloud that remained behind the cluster. Possibly 
this material was contained in a massive dense core within 
the cloud, which structure changed due to the process of 
cluster formation consuming the core. This scenario is con- 
sistent wit h the suggestion tha t NGC 6530 lies within an 
HII cavity (iMcCall et aLlfl990h . 

6. Discussion and Conclusions 

We developed a detailed study of the extinction in the di- 
rection of four clusters with ages < 5 Myr, located in star- 
forming regions. The aiming was to search for variable or 
anomalous extinction that possibly occurs for objects asso- 
ciated to dark clouds. In order to verify the characteristics 
of the clouds that coincide with the projected position of 
the cluster, we inspected the visual extinction maps and 
compared them to the far-IR emission detected in the di- 
rection of our sample. 

Different methods were used in order to improve previ- 
ous results on the extinction law Ry, providing means to 
identify the cause of anomalous extinction. In this compar- 
ative study we used the same photometric data base, ensur- 
ing similar observational conditions for the entire sample. 
Theoretical expressions for the extinction law were adopted 
to evaluate Ry based on the colour-excess E(B-V), which 
could be determined for some of the cluster members that 
have spectral type available in the literature. Variable ex- 
tinction law was clearly verified for NGC 2264, which pre- 
sented a bimodal distribution. Part of the cluster suffers 
normal extinction and another part shows large dispersion 
of Ry values, as illustrated in Fig. 3. However, these results 
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are not conclusive since few objects had spectral type well 
determined, particularly Berkeley 86. 

A more efficient analysis to estimate Ry and its de- 
pendence on wavelength was based on TCDs. The slope of 
the observed colours distribution, compared to the slope of 
ZAMS colours in Table 2, indicated anomalous extinction 
for NGC 6530, which gives Ry = 4.5 measured in the I 
band, and Rv — 6.2 in K band. The other clusters have 
normal extinction, in average. However, some members of 
NGC 2264 and NGC 2244 appear dispersed in the TCDs, 
indicating individual anomalous Ry, as illustrated in Figs. 
Al and A2. The TCD analysis confirms the results obtained 
from the colour-excess analysis, based on spectral type. 
Once there is no relation between the variation of Ry and 
the spatial distribution of the cluster members, we conclude 
that differences in the extinction law verified for some clus- 
ter members are not caused by environmental differences, 
but are probably due to circumstellar effect s. This is par- 
ticula rly in agreement with the results from lAlencar et ail 
(|2010f) that confirm the accretion activity in circumstellar 
discs of NGC 2264 members, for instance. 

In order to confirm if our results give reliable redden- 
ing correction, the Ry estimated from TCDs was used to 
determine the visual extinction that best fits the observed 
colours of the cluster members. A ZAMS fitting algorithm 
was developed to improve the search for Ay that minimizes 
the distance of the colours distribution to the reddened 
ZAMS curve, which was reproduced by Spline interpola- 
tion. The best fitting provided a mean value for E(B-V), 
which was compared to the results from other methods 
(Table 3). All the colour-excess estimations are in good 
agreement, excepting the E(B-V) derived from Ay maps 
mainly for NGC 2264 and NGC 6530. These clusters are 
not embedded on their respective clouds, which is consis- 
tent with the extinction levels found for the cluster mem- 
bers being lower than that measured on the extinction maps 
for the background stars. 

Fractal analysis was performed to investigate the sub- 
structures of the clouds and to compare them with statis- 
tical parameters of the cluster members distribution. The 
estimation of the Q parameter indicated that NGC 6530 
has a radially clustered spatial distribution, while the other 
clusters have fractally sub-clustered distribution of mem- 
bers. The fractal dimension (D2) measured on the Ay 
maps, which is related to the cloud geometry, was compared 
to Q measured on the clusters distributions. Excepting 
NGC 6530, the clusters have Q parameter compatible with 
the fractal dimension of the corresponding clouds, similar 
to the distribution of artificial data. We interpreted this 
result as a correlation between cloud structure and clus- 
ter members distribution, which is similar for Berkeley 86, 
NGC 2244 and NGC 2264. On the other hand, NGC 6530 
does not show this cloud-cluster relation, indicating that it 
was formed from a gas distribution more centrally concen- 
trated, different of the fractal sub-structures found in the 
remaining cloud. 

In fact, even under anomalous Ry — 4.5, NGC 6530 
does not suffer high extinction, as inferred from the average 
colour-excess (Ay ~ 1 mag), which is incompatible with the 
high levels of extinction shown in the dark clouds map or 
the far-IR emission maps. We conclude that anomalous ex- 
tinction in this case is not due to interstellar dense regions. 
A tentative explanation is the depletion of small grains due 
to evaporation that occurs under radiation from hot stars, 



which is consistent with the location of NGC 6530 in an HII 
cavity. Neither the circumstellar effects can be disregarded, 
since the grain growth is expected to occur in protoplan- 
etary discs. In this case, our results are also in agreement 
with the polarimetric results from McCall et al. (1990) that 
suggested the anomalous Ry is due to circumstellar effects. 

We also verified that the fractal analysis used to inves- 
tigate the cloud-cluster relation may give some indication 
about the scenario of the cluster formation. Once the role of 
this cloud-cluster relation is comparing the physical struc- 
ture of the cluster with its parental cloud, the interpretation 
of this relation may tell us about the original gas distribu- 
tion of the forming cloud and how the particular cloud may 
have evolved since the cluster was formed. It must be kept 
in mind that suggesting a cloud-cluster connection based 
on the correlation of Q with D2 is only speculative in the 
case of our sample, since more points are required to have 
a robust analysis. However, it is interesting to note the ob- 
served trend when comparing calculations performed for 
data sets from unrelated origins. This trend suggests a con- 
nection of stars clustering with sub-structures of the clouds, 
comparable to the artificial clouds and their derived clus- 
ters, studied by LWC11. For this reason, we conclude that 
NGC 6530 shows actual differences when compared with 
the other clusters, which is in agreement with previous re- 
sults from literature. 

An interesting perspective of this work is to extend our 
study to a larger number of clusters, mainly those that have 
not been studied as well as our sample has been. 
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Fig. A.l. TCDs for the clusters Berkeley 86 (left) and NGC 2244 (right). The solid line shows a normal extinction law, 
while the fitting for the cluster members is shown by a dotted line. Two extinction vectors (Ay = 2 mag) related to 
different values of Ry are plotted in the bottom panels, for illustration. Open circles indicate the members and crosses 
show the stars with P%< P m in- 
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Fig. A. 3. Top: Reddened ZAMS fitting to the cluster members with P%>P m ,„ (open circles). The position of the other 
stars of the sample is indicated by crosses. A full line shows the curve obtained with the fitting algorithm, while the 
dotted line shows the intrinsic U-B and B-V colours. Bottom: Curves representing the limit AE(B — V) = O.f 1 suggested 
by Burki (1975). From left to right: Berkeley 86, NGC 2244, NGC 2264, and NGC 6530. 
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